Alternative genome sequencing approaches of SARS-CoV-2 using Ion AmpliSeq Technology.

A thorough understanding of SARS-CoV-2 genetic features is compulsory to track the ongoing pandemic across multiple geographical locations of the world. Thermo Fisher Scientific USA has developed the Ion AmpliSeq SARS-CoV-2 Research Panel for the targeted sequencing of SARS-CoV-2 complete genome with high coverage and lower error rate. In this study an alternative approach of complete genome sequencing has been validated using different commercial sequencing kits to sequence the SARS-CoV-2. Amplification of cDNA with the SARS-CoV-2 primer pool was performed separately using two different master mixes: 2X environmental master mix (EM) and Platinum™ PCR SuperMix High Fidelity master mix (PM) instead of 5X Ion AmpliSeq™ HiFi Mix whereas NEBNext® Fast DNA Library Prep Set for Ion Torrent™ kit was used as an alternative to Ion AmpliSeq Library Kit Plus for other reagents. This study demonstrated a successful procedure to sequence the SARS-CoV-2 whole genome with average ∼2351 depth and 98.1% of total the reads aligned against the reference sequence (SARS-CoV-2, isolate Wuhan-Hu-1, complete genome). Although genome coverage varied, complete genomes were retrieved for both reagent sets with a reduced cost. This study proposed an alternative approach of high throughput sequencing using Ion torrent technology for the sequencing of SARS-CoV-2 in developing countries where sequencing facilities are low. This blended sequencing technique also offers a low cost protocol in developing countries like Bangladesh.

Novel insights into the co-selection of metal-driven antibiotic resistance in bacteria: a study of arsenic and antibiotic co-exposure.

The simultaneous development of antibiotic resistance in bacteria due to metal exposure poses a significant threat to the environment and human health. This study explored how exposure to both arsenic and antibiotics affects the ability of an arsenite oxidizer, Achromobacter xylosoxidans CAW4, to transform arsenite and its antibiotic resistance patterns. The bacterium was isolated from arsenic-contaminated groundwater in the Chandpur district of Bangladesh. We determined the minimum inhibitory concentration (MIC) of arsenite, cefotaxime, and tetracycline for A. xylosoxidans CAW4, demonstrating a multidrug resistance (MDR) trait. Following this determination, we aimed to mimic an environment where A. xylosoxidans CAW4 was exposed to both arsenite and antibiotics. We enabled the strain to grow in sub-MIC concentrations of 1 mM arsenite, 40 µg/mL cefotaxime, and 20 µg/mL tetracycline. The expression dynamics of the arsenite oxidase (aioA) gene in the presence or absence of antibiotics were analyzed. The findings indicated that simultaneous exposure to arsenite and antibiotics adversely affected the bacteria's capacity to metabolize arsenic. However, when arsenite was present in antibiotics-containing media, it promoted bacterial growth. The study observed a global downregulation of the aioA gene in arsenic-antibiotic conditions, indicating the possibility of increased susceptibility through co-resistance across the entire bacterial population of the environment. This study interprets that bacterial arsenic-metabolizing ability can rescue the bacteria from antibiotic stress, further disseminating environmental cross-resistance. Therefore, the co-selection of metal-driven antibiotic resistance in bacteria highlights the need for effective measures to address this emerging threat to human health and the environment.

Unveiling distinct genetic features in multidrug-resistant Escherichia coli isolated from mammary tissue and gut of mastitis induced mice.

Escherichia coli is one of the major pathogens causing mastitis in lactating mammals. We hypothesized that E. coli from the gut and mammary glands may have similar genomic characteristics in the causation of mastitis. To test this hypothesis, we used whole genome sequencing to analyze two multidrug resistant E. coli strains isolated from mammary tissue (G2M6U) and fecal sample (G6M1F) of experimentally induced mastitis mice. Both strains showed resistance to multiple (>7) antibiotics such as oxacillin, aztreonam, nalidixic acid, streptomycin, gentamicin, cefoxitin, ampicillin, tetracycline, azithromycin and nitrofurantoin. The genome of E. coli G2M6U had 59 antimicrobial resistance genes (ARGs) and 159 virulence factor genes (VFGs), while the E. coli G6M1F genome possessed 77 ARGs and 178 VFGs. Both strains were found to be genetically related to many E. coli strains causing mastitis and enteric diseases originating from different hosts and regions. The G6M1F had several unique ARGs (e.g., QnrS1, sul2, tetA, tetR, emrK, blaTEM-1/105, and aph(6)-Id, aph(3″)-Ib) conferring resistance to certain antibiotics, whereas G2M6U had a unique heat-stable enterotoxin gene (astA) and 7192 single nucleotide polymorphisms. Furthermore, there were 43 and 111 unique genes identified in G2M6U and G6M1F genomes, respectively. These results indicate distinct differences in the genomic characteristics of E. coli strain G2M6U and G6M1F that might have important implications in the pathophysiology of mammalian mastitis, and treatment strategies for mastitis in dairy animals.

Genomic features and pathophysiological impact of a multidrug-resistant Staphylococcus warneri variant in murine mastitis.

Non-aureus staphylococci (NAS) represent a major etiological agent in dairy animal mastitis, yet their role and impact remain insufficiently studied. This study aimed to elucidate the genomic characteristics of a newly identified multidrug-resistant NAS strain, specifically Staphylococcus warneri G1M1F, isolated from murine feces in an experimental mastitis model. Surprisingly, NAS species accounted for 54.35 % of murine mastitis cases, with S. warneri being the most prevalent at 40.0 %. S. warneri G1M1F exhibited resistance to 10 major antibiotics. Whole-genome sequencing established a genetic connection between G1M1F and S. warneri strains isolated previously from various sources including mastitis milk in dairy animals, human feces and blood across diverse geographical regions. Genomic analysis of S. warneri G1M1F unveiled 34 antimicrobial resistance genes (ARGs), 30 virulence factor genes (VFGs), and 278 metabolic features. A significant portion of identified ARGs (64 %) conferred resistance through antibiotic efflux pumps, while VFGs primarily related to bacterial adherence and biofilm formation. Inoculation with G1M1F in mice resulted in pronounced inflammatory lesions in mammary and colon tissues, indicating pathogenic potential. Our findings highlight distinctive genomic traits in S. warneri G1M1F, signifying the emergence of a novel multidrug-resistant NAS variant. These insights contribute to understanding NAS-related mastitis pathophysiology and inform strategies for effective treatment in dairy animals.

Metagenomic and culture-dependent approaches unveil active microbial community and novel functional genes involved in arsenic mobilization and detoxification in groundwater.

BACKGROUND: Arsenic (As) and its species are major pollutants in ecological bodied including groundwater in Bangladesh rendering serious public health concern. Bacteria with arsenotrophic genes have been found in the aquifer, converting toxic arsenite [As (III)] to less toxic arsenate [As (V)] that is easily removed using chemical and biological trappers. In this study, genomic and metagenomic approaches parallel to culture-based assay (Graphical abstract) have made it possible to decipher phylogenetic diversity of groundwater arsenotrophic microbiomes along with elucidation of their genetic determinants. RESULTS: Seventy-two isolates were retrieved from six As-contaminated (average As concentration of 0.23 mg/L) groundwater samples from Munshiganj and Chandpur districts of Bangladesh. Twenty-three isolates harbored arsenite efflux pump (arsB) gene with high abundance, and ten isolates possessing arsenite oxidase (aioA) gene, with a wide range of minimum inhibitory concentration, MICAs (2 to 32 mM), confirming their role in arsenite metabolism. There was considerable heterogeneity in species richness and microbial community structure. Microbial taxa from Proteobacteria, Firmicutes and Acidobacteria dominated these diversities. Through these combinatorial approaches, we have identified potential candidates such as, Pseudomonas, Acinetobacter, Stenotrophomonas, Achromobacter, Paraburkholderia, Comamonas and Klebsiella and associated functional genes (arsB, acr3, arsD, arsH, arsR) that could significantly contribute to arsenite detoxification, accumulation, and immobilization. CONCLUSIONS: Culture-dependent and -independent shotgun metagenomic investigation elucidated arsenotrophic microbiomes and their functions in As biogeochemical transformation. These findings laid a foundation for further large-scale researches on the arsenotrophic microbiomes and their concurrent functions in As biogeochemical transformation in As-contaminated areas of Bangladesh and beyond.

Emergence of a novel sublineage, MYMBD21 under SA-2018 lineage of Foot-and-Mouth Disease Virus serotype O in Bangladesh.

Foot-and-Mouth Disease (FMD) hinders the growth of the livestock industry in endemic countries like Bangladesh. The management and prevention of FMD are severely impacted by the high mutation rate and subsequent frequent generation of newer genotypes of the causative agent, Foot-and-Mouth Disease Virus (FMDV). The current study was conducted in nine districts of Bangladesh during 2019-21 to characterize the circulating FMDV strains based on the VP1 sequence analysis, the major antigenic recognition site providing serotype specificity and high variability of FMDV. This study detected the first emergence of the SA-2018 lineage in Bangladesh along with the predominance of Ind-2001e (or Ind-2001BD1) sublineage of ME-SA topotype under serotype O during 2019-21. The mutational spectrum, evolutionary divergence analysis and multidimensional plotting confirmed the isolates collected from Mymensingh districts, designated as MYMBD21 as a novel sublineage under the SA-2018 lineage. Analysis of the amino acid sequence revealed several changes in the G-H loop, B-C loop and C-terminal region of VP1, revealing a 12-13% divergence from the existing vaccine strains and a 95% VP1 protein homology, with most of the mutations potentially considerable as vaccine escape mutations, evidenced by three-dimensional structural analysis. This is the first report on the emergence of the SA-2018 lineage of ME-SA topotype of FMDV serotype O in Bangladesh, as well as a possible mutational trend towards the emergence of a distinct sublineage under SA-2018 lineage, which calls for in-depth genome-wide analysis and monitoring of the FMD situation in the country to implement a strategic vaccination and effective FMD control program.

Genomic variability correlates with biofilm phenotypes in multidrug resistant clinical isolates of Pseudomonas aeruginosa.

The multifactorial nature of Pseudomonas aeruginosa biofilm development and genomic variabilities implicates its resistance to conventional antimicrobials and virulence. Therefore, genetic determinants need to be extensively studied to block the early steps of biofilm or already formed biofilms. In this study, a total of 20 multidrug resistant (MDR) clinical P. aeruginosa isolates were evaluated for their biofilm forming abilities and related genes. Of the isolates tested, all of them showed surface attachment tendencies in nutrient limiting conditions, and classified as strong (SBF = 45%), moderate (MBF = 30%) and weak (WBF = 25%) biofilm formers. Complete genome sequencing of representative strong (DMC-27b), moderate (DMC-20c) and weak biofilm former (DMC-30b) isolates was performed. Analysis of biofilm related genes in the sequenced genomes revealed that, 80 of the 88 biofilm related genes possess 98-100% sequence identity to the reference PAO1 strain. Complete and partial sequence data of LecB proteins from tested isolates indicate that isolates containing PA14-like LecB sequences produced strong biofilms. All of the 7 pel operon protein coding genes in weak biofilm former isolate 30b showed significant nucleotide sequence variation with other tested isolates, and their corresponding proteins are 99% identical with the pel operon proteins of PA7. Bioinformatics analyses identified divergent sequence and structural features that separate PA7 like pel operon proteins from reference PAO1-like pel operon. Congo red and pellicle forming assays revealed that the sequence and structure variations may have interfered with the Pel production pathway and resulted in impaired Pel production in isolate 30b that has a PA7 like pel operon. Expression analysis also showed that both pelB and lecB genes were about 5 to 6 folds upregulated after 24 h in SBF 27b in comparison with WBF 30b. Our findings indicate significant genomic divergence in biofilm related genes of P. aeruginosa strains that affect their biofilm phenotypes.

Temporal dynamics and fatality of SARS-CoV-2 variants in Bangladesh.

Background and Aims: Since the beginning of the SARS-CoV-2 pandemic, multiple new variants have emerged posing an increased risk to global public health. This study aimed to investigate SARS-CoV-2 variants, their temporal dynamics, infection rate (IFR) and case fatality rate (CFR) in Bangladesh by analyzing the published genomes. Methods: We retrieved 6610 complete whole genome sequences of the SARS-CoV-2 from the GISAID (Global Initiative on Sharing all Influenza Data) platform from March 2020 to October 2022, and performed different in-silico bioinformatics analyses. The clade and Pango lineages were assigned by using Nextclade v2.8.1. SARS-CoV-2 infections and fatality data were collected from the Institute of Epidemiology Disease Control and Research (IEDCR), Bangladesh. The average IFR was calculated from the monthly COVID-19 cases and population size while average CFR was calculated from the number of monthly deaths and number of confirmed COVID-19 cases. Results: SARS-CoV-2 first emerged in Bangladesh on March 3, 2020 and created three pandemic waves so far. The phylogenetic analysis revealed multiple introductions of SARS-CoV-2 variant(s) into Bangladesh with at least 22 Nextstrain clades and 107 Pangolin lineages with respect to the SARS-CoV-2 reference genome of Wuhan/Hu-1/2019. The Delta variant was detected as the most predominant (48.06%) variant followed by Omicron (27.88%), Beta (7.65%), Alpha (1.56%), Eta (0.33%) and Gamma (0.03%) variant. The overall IFR and CFR from circulating variants were 13.59% and 1.45%, respectively. A time-dependent monthly analysis showed significant variations in the IFR (p = 0.012, Kruskal-Wallis test) and CFR (p = 0.032, Kruskal-Wallis test) throughout the study period. We found the highest IFR (14.35%) in 2020 while Delta (20A) and Beta (20H) variants were circulating in Bangladesh. Remarkably, the highest CFR (1.91%) from SARS-CoV-2 variants was recorded in 2021. Conclusion: Our findings highlight the importance of genomic surveillance for careful monitoring of variants of concern emergence to interpret correctly their relative IFR and CFR, and thus, for implementation of strengthened public health and social measures to control the spread of the virus. Furthermore, the results of the present study may provide important context for sequence-based inference in SARS-CoV-2 variant(s) evolution and clinical epidemiology beyond Bangladesh.

Evaluation of RT-PCR assays for detection of SARS-CoV-2 variants of concern.

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) pandemic has been considered with great importance on correct screening procedure. The detection efficiency of recent variants of concern were observed by comparing 5 commercial RT-PCR kits and a SYBR-green method developed and validated in our laboratory. The RNA was extracted from nasopharyngeal samples from suspected COVID-19 patients and RT-PCR assay was performed according to the instruction of the respective manufacturers. The specificity and sensitivity of Maccura kit was 81.8% and 82.5%, A*Star kit was 100% and 75.4%, Da An Gene kit was 100% and 68.4%, Sansure kit was 54.5% and 91.2% and TaqPath kit was 100% and 70.2% respectively. Our in house SYBR-Green method showed a consistent detection result with 90.9% specificity and 91.2% sensitivity. We also found that detection kits targeting more genes showed better accuracy which facilitates less false positive results (< 20%). Our study found a significant difference (p < 0.005) in Ct value reported for common target genes shared by the RT-PCR kits in relation with different variants of COVID-19 infection. Recent variants of concerns contain more than 30 mutations in the spike proteins including 2 deletion and a unique insertion mutation by which makes detection of these variants difficult and these facilitates the variants to escape from being detected.

Transcriptome analysis reveals increased abundance and diversity of opportunistic fungal pathogens in nasopharyngeal tract of COVID-19 patients.

We previously reported that SARS-CoV-2 infection reduces human nasopharyngeal commensal microbiomes (bacteria, archaea and commensal respiratory viruses) with inclusion of pathobionts. This study aimed to assess the possible changes in the abundance and diversity of resident mycobiome in the nasopharyngeal tract (NT) of humans due to SARS-CoV-2 infections. Twenty-two (n = 22) nasopharyngeal swab samples (including COVID-19 = 8, Recovered = 7, and Healthy = 7) were collected for RNA-sequencing followed by taxonomic profiling of mycobiome. Our analyses indicate that SARS-CoV-2 infection significantly increased (p < 0.05, Wilcoxon test) the population and diversity of fungi in the NT with inclusion of a high proportion of opportunistic pathogens. We detected 863 fungal species including 533, 445, and 188 species in COVID-19, Recovered, and Healthy individuals, respectively that indicate a distinct mycobiome dysbiosis due to the SARS-CoV-2 infection. Remarkably, 37% of the fungal species were exclusively associated with SARS-CoV-2 infection, where S. cerevisiae (88.62%) and Phaffia rhodozyma (10.30%) were two top abundant species. Likewise, Recovered humans NT samples were predominated by Aspergillus penicillioides (36.64%), A. keveii (23.36%), A. oryzae (10.05%) and A. pseudoglaucus (4.42%). Conversely, Nannochloropsis oceanica (47.93%), Saccharomyces pastorianus (34.42%), and S. cerevisiae (2.80%) were the top abundant fungal species in Healthy controls nasal swabs. Importantly, 16% commensal fungal species found in the Healthy controls were not detected in either COVID-19 patients or when they were cured from COVID-19 (Recovered). We also detected several altered metabolic pathways correlated with the dysbiosis of fungal mycobiota in COVID-19 patients. Our results suggest that SARS-CoV-2 infection causes significant dysbiosis of mycobiome and related metabolic functions possibly play a determining role in the progression of SARS-CoV-2 pathogenesis. These findings might be helpful for developing mycobiome-based diagnostics, and also devising appropriate therapeutic regimens including antifungal drugs for prevention and control of concurrent fungal coinfections in COVID-19 patients.

Genomic diversity and molecular epidemiology of a multidrug-resistant Pseudomonas aeruginosa DMC30b isolated from a hospitalized burn patient in Bangladesh.

OBJECTIVES: Pseudomonas aeruginosa is a key opportunistic pathogen causing a wide range of community- and hospital-acquired infections in immunocompromised or catheterized patients. Here, we report the complete genome sequence of a multidrug-resistant (MDR) P. aeruginosa DMC30b to elucidate the genetic diversity, molecular epidemiology, and underlying mechanisms for antimicrobial resistance and virulence. METHODS: P. aeruginosa DMC30b was isolated from septic wound swab of a severe burn patient. Whole-genome sequencing was performed under Ion Torrent platform. The genome was assembled using the SPAdes v. 3.12.01 in an integrated Genome Analysis Platform for Ion Torrent sequence data. The genome was annotated using the National Center for Biotechnology Information (NCBI) Prokaryotic Genome Annotation Pipeline. In-silico predictions of antimicrobial resistance genes, virulence factor genes, and metabolic functional potentials were performed using different curated bioinformatics tools. RESULTS: P. aeruginosa DMC30b was found as a MDR strain and belonged to sequence type 244 (ST244). The complete genome size is 6 994 756 bp with a coverage of 76.76x, guanine-cytosine content of 65.7% and a Benchmarking Universal Single-Copy Orthologs score of 100. The genome of P. aeruginosa DMC30b harboured two predicted plasmid replicons (e,g. IncP-6; 78 007 bp and ColRNAI; 9359 bp), 35 resistomes (antimicrobial resistance genes) conferring resistance to 18 different antibiotics (including four beta-lactam classes), and 214 virulence factor genes. It was identified as the 167th ST244 strain among ∼ 5800 whole-genome sequences of P. aeruginosa available in the NCBI database. CONCLUSION: The MDR P. aeruginosa DMC30b was identified as the 167th ST244 complete genome to be submitted to the NCBI, and the first ST244 isolate sequenced from Bangladesh. The complete genome data with high genetic diversity and underlying mechanisms for antimicrobial resistance and virulence of P. aeruginosa DMC30b will aid in understanding the evolution and phylogeny of such high-risk clones and provide a solid basis for further research on MDR or extensively drug resistant strains.

Induction of mastitis by cow-to-mouse fecal and milk microbiota transplantation causes microbiome dysbiosis and genomic functional perturbation in mice.

BACKGROUND: Mastitis pathogenesis involves a wide range of opportunistic and apparently resident microorganims including bacteria, viruses and archaea. In dairy animals, microbes reside in the host, interact with environment and evade the host immune system, providing a potential for host-tropism to favor mastitis pathogenesis. To understand the host-tropism phenomena of bovine-tropic mastitis microbiomes, we developed a cow-to-mouse mastitis model. METHODS: A cow-to-mouse mastitis model was established by fecal microbiota transplantation (FMT) and milk microbiota transplantation (MMT) to pregnant mice to assess microbiome dysbiosis and genomic functional perturbations through shotgun whole metagenome sequencing (WMS) along with histopathological changes in mice mammary gland and colon tissues. RESULTS: The cow-to-mouse FMT and MMT from clinical mastitis (CM) cows induced mastitis syndromes in mice as evidenced by histopathological changes in mammary gland and colon tissues. The WMS of 24 samples including six milk (CM = 3, healthy; H = 3), six fecal (CM = 4, H = 2) samples from cows, and six fecal (CM = 4, H = 2) and six mammary tissue (CM = 3, H = 3) samples from mice generating 517.14 million reads (average: 21.55 million reads/sample) mapped to 2191 bacterial, 94 viral and 54 archaeal genomes. The Kruskal-Wallis test revealed significant differences (p = 0.009) in diversity, composition, and relative abundances in microbiomes between CM- and H-metagenomes. These differences in microbiome composition were mostly represented by Pseudomonas aeruginosa, Lactobacillus crispatus, Klebsiella oxytoca, Enterococcus faecalis, Pantoea dispersa in CM-cows (feces and milk), and Muribaculum spp., Duncaniella spp., Muribaculum intestinale, Bifidobacterium animalis, Escherichia coli, Staphylococcus aureus, Massilia oculi, Ralstonia pickettii in CM-mice (feces and mammary tissues). Different species of Clostridia, Bacteroida, Actinobacteria, Flavobacteriia and Betaproteobacteria had a strong co-occurrence and positive correlation as the indicator species of murine mastitis. However, both CM cows and mice shared few mastitis-associated microbial taxa (1.14%) and functional pathways regardless of conservation of mastitis syndromes, indicating the higher discrepancy in mastitis-associated microbiomes among lactating mammals. CONCLUSIONS: We successfully induced mastitis by FMT and MMT that resulted in microbiome dysbiosis and genomic functional perturbations in mice. This study induced mastitis in a mouse model through FMT and MMT, which might be useful for further studies- focused on pathogen(s) involved in mastitis, their cross-talk among themselves and the host.

Development and validation of cost-effective one-step multiplex RT-PCR assay for detecting the SARS-CoV-2 infection using SYBR Green melting curve analysis.

TaqMan probe-based commercial real-time (RT) PCR kits are expensive but most frequently used in COVID-19 diagnosis. The unprecedented scale of SARS-CoV-2 infections needs to meet the challenge of testing more persons at a reasonable cost. This study developed a simple and cost-effective alternative diagnostic method based on melting curve analysis of SYBR green multiplex assay targeting two virus-specific genes along with a host-specific internal control. A total of 180 randomly selected samples portioning into two subsets based on crude and high-quality RNA extraction were used to compare this assay with a nationwide available commercial kit (Sansure Biotech Inc., (Hunan, China)), so that we could analyze the variation and validity of this in-house developed method. Our customized-designed primers can specifically detect the viral RNA likewise Sansure. We separately optimized SYBR Green RT-PCR reaction of N, E, S, and RdRp genes based on singleplex melting curve analysis at the initial stage. After several rounds of optimization on multiplex assays of different primer combinations, the optimized method finally targeted N and E genes of the SARS-CoV-2 virus, together with the ß-actin gene of the host as an internal control. Comparing with the Sansure commercial kit, our proposed assay provided up to 97% specificity and 93% sensitivity. The cost of each sample processing ranged between ~2 and ~6 USD depending on the purification level of extracted RNA template. Overall, this one-step and one-tube method can revolutionize the COVID-19 diagnosis in low-income countries.

Circulating Phylotypes of White Spot Syndrome Virus in Bangladesh and Their Virulence.

White Spot Syndrome Virus (WSSV) has emerged as one of the most prevalent and lethal viruses globally and infects both shrimps and crabs in the aquatic environment. This study aimed to investigate the occurrence of WSSV in different ghers of Bangladesh and the virulence of the circulating phylotypes. We collected 360 shrimp (Penaeus monodon) and 120 crab (Scylla sp.) samples from the south-east (Cox's Bazar) and south-west (Satkhira) coastal regions of Bangladesh. The VP28 gene-specific PCR assays and sequencing revealed statistically significant (p < 0.05, Kruskal-Wallis test) differences in the prevalence of WSSV in shrimps and crabs between the study areas (Cox's Bazar and Satkhira) and over the study periods (2017-2019). The mean Log load of WSSV varied from 8.40 (Cox's Bazar) to 10.48 (Satkhira) per gram of tissue. The mean values for salinity, dissolved oxygen, temperature and pH were 14.71 ± 0.76 ppt, 3.7 ± 0.1 ppm, 34.11 ± 0.38 °C and 8.23 ± 0.38, respectively, in the WSSV-positive ghers. The VP28 gene-based phylogenetic analysis showed an amino-acid substitution (E→G) at the 167th position in the isolates from Cox's Bazar (referred to as phylotype BD2) compared to the globally circulating one (BD1). Shrimp PL artificially challenged with BD1 and BD2 phylotypes with filtrates of tissue containing 0.423 × 109 copies of WSSV per mL resulted in a median LT50 value of 73 h and 75 h, respectively. The in vivo trial showed higher mean Log WSSV copies (6.47 ± 2.07 per mg tissue) in BD1-challenged shrimp PL compared to BD2 (4.75 ± 0.35 per mg tissue). Crabs infected with BD1 and BD2 showed 100% mortality within 48 h and 62 h of challenge, respectively, with mean Log WSSV copies of 12.06 ± 0.48 and 9.95 ± 0.37 per gram tissue, respectively. Moreover, shrimp antimicrobial peptides (AMPs), penaeidin and lysozyme expression were lower in the BD1-challenged group compared to BD2 challenged shrimps. These results collectively demonstrated that relative virulence properties of WSSV based on mortality rate, viral load and expression of host immune genes in artificially infected shrimp PL could be affected by single aa substitution in VP28.

Dominant clade-featured SARS-CoV-2 co-occurring mutations reveal plausible epistasis: An in silico based hypothetical model.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved into eight fundamental clades with four of these clades (G, GH, GR, and GV) globally prevalent in 2020. To explain plausible epistatic effects of the signature co-occurring mutations of these circulating clades on viral replication and transmission fitness, we proposed a hypothetical model using in silico approach. Molecular docking and dynamics analyses showed the higher infectiousness of a spike mutant through more favorable binding of G614 with the elastase-2. RdRp mutation p.P323L significantly increased genome-wide mutations (p < 0.0001), allowing for more flexible RdRp (mutated)-NSP8 interaction that may accelerate replication. Superior RNA stability and structural variation at NSP3:C241T might impact protein, RNA interactions, or both. Another silent 5'-UTR:C241T mutation might affect translational efficiency and viral packaging. These four G-clade-featured co-occurring mutations might increase viral replication. Sentinel GH-clade ORF3a:p.Q57H variants constricted the ion-channel through intertransmembrane-domain interaction of cysteine(C81)-histidine(H57). The GR-clade N:p.RG203-204KR would stabilize RNA interaction by a more flexible and hypo-phosphorylated SR-rich region. GV-clade viruses seemingly gained the evolutionary advantage of the confounding factors; nevertheless, N:p.A220V might modulate RNA binding with no phenotypic effect. Our hypothetical model needs further retrospective and prospective studies to understand detailed molecular events and their relationship to the fitness of SARS-CoV-2.

Genomic analysis of SARS-CoV-2 variants of concern identified from the ChAdOx1 nCoV-19 immunized patients from Southwest part of Bangladesh.

BACKGROUND: Bangladesh introduced ChAdOx1 nCoV-19 since February, 2021 and in six months, only a small population (12.8%) received either one or two dose of vaccination like other low-income countries. The COVID-19 infections were continued to roll all over the places although the information on genomic variations of SARS-CoV-2 between both immunized and unimmunized group was unavailable. The objective of this study was to compare the proportion of immune escaping variants between those groups. METHODS: A total of 4718 nasopharygeal samples were collected from March 1 until April 15, 2021, of which, 834 (18%) were SARS-CoV-2 positive. The minimum sample size was calculated as 108 who were randomly selected for telephone interview and provided consent. The prevalence of SARS-CoV-2 variants and disease severity among both immunized and unimmunized groups was measured. A total of 63 spike protein sequences and 14 whole-genome sequences were performed from both groups and phylogenetic reconstruction and mutation analysis were compared. RESULTS: A total of 40 respondents (37%, N = 108) received single-dose and 2 (2%) received both doses of ChAdOx1 nCoV-19 vaccine, which significantly reduce dry cough, loss of appetite and difficulties in breathing compared to none. There was no significant difference in hospitalization, duration of hospitalization or reduction of other symptoms like running nose, muscle pain, shortness of breathing or generalized weakness between immunized and unimmunized groups. Spike protein sequence assumed 21 (87.5%) B.1.351, one B.1.526 and two 20B variants in immunized group compared to 27 (69%) B.1.351, 5 (13%) B.1.1.7, 4 (10%) 20B, 2 B.1.526 and one B.1.427 variant in unimmunized group. Whole genome sequence analysis of 14 cases identified seven B.1.351 Beta V2, three B.1.1.7 Alpha V1, one B.1.526 Eta and the rest three 20B variants. CONCLUSION: Our study observed that ChAdOx1 could not prevent the new infection or severe COVID-19 disease outcome with single dose while the infections were mostly caused by B.1.351 variants in Bangladesh.

First report from Bangladesh on genetic diversity of multidrug-resistant Pasteurella multocida type B:2 in fowl cholera.

BACKGROUND AND AIM: Fowl cholera (FC) caused by Pasteurella multocida is a highly contagious bacterial disease of global importance for poultry production. The severity and incidence of FC caused by P. multocida may vary considerably depending on several factors associated with the host (including species and age of infected birds), the environment, and the bacterial strain. This study aimed to investigate the genetic diversity of multidrug-resistant P. multocida strains isolated from FC outbreaks in laying hens from commercial farms of Bangladesh. MATERIALS AND METHODS: We collected 57 samples of suspected FC, including 36 live and 21 dead laying hens. P. multocida isolates were characterized by biochemical and molecular-biological methods. RESULTS: Twenty-two strains of P. multocida were isolated from these samples through phenotypic and genotypic characterization. The strains were grouped into two distinct random amplification of polymorphic DNA (RAPD) biotypes harboring a range of pathogenic genes; exbB, ompH, ptfA, nanB, sodC, and hgbA. In this study, 90.90% and 81.82% P. multocida strains were multidrug-resistant and biofilm formers, respectively. Whole-genome sequencing of the two representative RAPD phylotypes confirmed as P. multocida type B: L2:ST122, harboring a number of virulence factors-associated genes (VFGs), and antimicrobial resistance (AMR) genes (ARGs). In addition, pan-genome analysis revealed 90 unique genes in the genomes of P. multocida predicted to be associated with versatile metabolic functions, pathogenicity, virulence, and AMR. CONCLUSION: This is first-ever report on the association of P. multocida genotype B: L2:ST122 and related VFGs and ARGs in the pathogenesis of FC in laying hens. This study also provides a genetic context for future researches on the evolutionary diversity of P. multocida strains and their host adaptation.

Draft genome sequence of a carbapenem-resistant clinical Acinetobacter baumannii revealing co-existence of four classes of ß-lactamases.

OBJECTIVES: Nosocomial carbapenem-resistant Acinetobacter baumannii is a challenge in the treatment of intensive care unit (ICU) patients. The presence of antimicrobial resistance genes (ARGs) and mobile genetic elements can further complicate effects to combat antibiotic resistance in this post-antibiotic era. The aim of this study was to analyse the molecular basis of carbapenem resistance in A. baumannii strain DMC-32a. METHODS: Strain DMC-32a, isolated from a wound swab of an ICU patient, was screened phenotypically and genotypically for carbapenem resistance. The isolate was subjected whole-genome sequencing (WGS) to investigate the resistance mechanisms. RESULTS: Strain DMC-32a was resistant to all tested antibiotics belonging to seven classes, except for the polymyxins. MICs determined against imipenem and meropenem were 512 mg/L and >512 mg/L, respectively. Gene-specific PCR confirmed the presence of blaNDM-1 and intI1. WGS confirmed the isolate as sequence type ST1053, the presence of four classes of ß-lactamases [A (blaPER-7, blaTEM-116), B (blaNDM-1), C (blaADC-26) and D (blaOXA-23, blaOXA-51)], and four genes encoding aminoglycoside-modifying enzymes [aph(6)-Id, aph(3'')-Ib, aph(3')-VI and ant(3'')-IIc] conferring resistance to streptomycin, amikacin, kanamycin and spectinomycin. The analysis also confirmed the presence of a class 1 integron gene cassette harbouring ARGs to rifamycin (arr-2) and phenicols (cmlA5). No plasmid replicon was found from the sequence data. CONCLUSION: The co-existence of four ß-lactamase classes in A. baumannii DMC-32a has not been reported previously from Bangladesh and not in this ST elsewhere. The emergence of such a nosocomial pathogen creates the need for effective therapeutics for critically-ill patients and for controlling hospital-acquired infections.

Microbiome signature and diversity regulates the level of energy production under anaerobic condition.

The microbiome of the anaerobic digester (AD) regulates the level of energy production. To assess the microbiome diversity and composition in different stages of anaerobic digestion, we collected 16 samples from the AD of cow dung (CD) origin. The samples were categorized into four groups (Group-I, Group-II, Group-III and Group-IV) based on the level of energy production (CH4%), and sequenced through whole metagenome sequencing (WMS). Group-I (n = 2) belonged to initial time of energy production whereas Group-II (n = 5), Group-III (n = 5), and Group-IV (n = 4) had 21-34%, 47-58% and 71-74% of CH4, respectively. The physicochemical analysis revealed that level of energy production (CH4%) had significant positive correlation with digester pH (r = 0.92, p < 0.001), O2 level (%) (r = 0.54, p < 0.05), and environmental temperature (°C) (r = 0.57, p < 0.05). The WMS data mapped to 2800 distinct bacterial, archaeal and viral genomes through PathoScope (PS) and MG-RAST (MR) analyses. We detected 768, 1421, 1819 and 1774 bacterial strains in Group-I, Group-II, Group-III and Group-IV, respectively through PS analysis which were represented by Firmicutes, Bacteroidetes, Proteobacteria, Actinobacteria, Spirochaetes and Fibrobacteres phyla (> 93.0% of the total abundances). Simultaneously, 343 archaeal strains were detected, of which 95.90% strains shared across four metagenomes. We identified 43 dominant species including 31 bacterial and 12 archaeal species in AD microbiomes, of which only archaea showed positive correlation with digester pH, CH4 concentration, pressure and temperature (Spearman correlation; r > 0.6, p < 0.01). The indicator species analysis showed that the species Methanosarcina vacuolate, Dehalococcoides mccartyi, Methanosarcina sp. Kolksee and Methanosarcina barkeri were highly specific for energy production. The correlation network analysis showed that different strains of Euryarcheota and Firmicutes phyla exhibited significant correlation (p = 0.021, Kruskal-Wallis test; with a cutoff of 1.0) with the highest level (74.1%) of energy production (Group-IV). In addition, top CH4 producing microbiomes showed increased genomic functional activities related to one carbon and biotin metabolism, oxidative stress, proteolytic pathways, membrane-type-1-matrix-metalloproteinase (MT1-MMP) pericellular network, acetyl-CoA production, motility and chemotaxis. Importantly, the physicochemical properties of the AD including pH, CH4 concentration (%), pressure, temperature and environmental temperature were found to be positively correlated with these genomic functional potentials and distribution of ARGs and metal resistance pathways (Spearman correlation; r > 0.5, p < 0.01). This study reveals distinct changes in composition and diversity of the AD microbiomes including different indicator species, and their genomic features that are highly specific for energy production.